Topological Disorder in Spin Models on Hierarchical Lattices. 
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A general approach for the description of spin systems on hierarchial lattices with coordination 
number g as a dynamical variable is proposed. The ferromagnetic Ising model on the Bethe lattice 
was studied as a simple example demonstrating our method. The annealed and partly annealed 
versions of disorder concerned with the lattice coordination number are invented and discussed. Re- 
current relations are obtained for the evaluation of magnetization. The magnetization is calculated 
for the particular disorder choices, q — 2,3 and q = 3,4. A nontrivial localization of critical point is 
revealed. 
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I. INTRODUCTION 

The hierarchical models in statistical physics can be 
used to exactly describe the critical phenomena. For 
these purposes only approximate or numerical methods 
exist in most statistical mechanical models. Therefore 
complicated models representing nontrivial critical be- 
havior and demonstrating analytical solution may be rec- 
ognized as very important. A good example of such sys- 
tems are spin models on hierarchical lattices. Besides the 
general interest related to the critical phenomena, spin 
models on hierarchical lattices are also formulated for 
description of the concrete physical systems. The inter- 
esting applications are Bethe approximation in polymer 
physics, and many others [H-Q. 

The main purpose of this work is a study of spin models 
on hierarchical lattices with so called topological disor- 
der. The problem of disorder for the hierarchical models 
is known and has been discussed in series of works p^ - 
[T3| . However, disorders concerned with couplings and 
magnetic fields are considered mainly. So called topo- 
logical disorder, describing by the changing lattice coor- 
dination number, is investigated less. In our definition 
topological disorder represents non-homogeneity of the 
lattice structure and may be even non-connected to the 
spins dynamics. A more precise definition of the topo- 
logical disorder in context of the hierarchical models is 
given below within a concrete example. 

Although all results in our article are related to the 
physics and have direct meaning in a physical sense, it 
is interesting to discuss their possible applications in the 
context of social dynamics. As is known, in the recent 
decades there is a rapidly growing interest to the appli- 
cation of the statistical physics to interdisciplinary fields 
(as biology, computer science, social dynamics and e.t.c). 
Considering the social dynamics we deal not with phys- 
ical particles but with humans or crowds. The human 
behavior is not rather well understood and depends on 
many hidden factors. Nevertheless, one may suppose 
that in the statistical limit the new statistical laws would 
be obtained, which are less dependent on the individuals. 
As it was shown in recent decade, this conjecture is cor- 



rect. The transitions from disorder to order are discov- 
ered describing the social dynamics from the statistical 
physics point of view. We may bring as the examples 
of those transitions spontaneous formation of a common 
language, culture or the emergence of consensus about a 
specific issue • There are also examples of scaling and 
universality. 

As an illustration of our approach to spin systems on 
hierarchical lattices with changing coordination number 
q we consider ferromagnetic Ising model on the Bethe 
lattice as a simplest possible example. 

The article is organized as follows. In section |lT] a sim- 
ple ferromagnetic Ising model with topological disorder 
on hierarchical lattice is formulated. A more precise and 
clear definition of the topological disorder in context of 
the hierarchical lattices is given. The connection between 
proposed model and description of social dynamics is dis- 
cussed in detail. 

In section Hill annealed and partly annealed topological 
disorders are defined. A general method for those mathe- 
matical description is proposed. The annealed version of 
the topological disorder is studied for the particular real- 
izations of the lattice structure. The critical temperature 
is found and analyzed. 

Partly annealed and quenched topological disorders are 
described within ferromagnetic Ising model in section lTVl 
The same particular realizations of the topological disor- 
der as for the annealed one are discussed. The critical 
point is obtained. 

Main observations are presented in conclusion. 



II. ISING MODEL ON BETHE LATTICE WITH 
TOPOLOGICAL DISORDER 

The one of the simplest recursive lattices is a Cayley 
tree. Constructing Cayley tree one starts with a central 
site O. After that we connect the central site O by links 
with q new other sites. This procedure of q new sites 
connection by links to the each site originated on the 
previous step repeats n times. Following this procedure 
one gets the recursive connected graph, containing no cy- 
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cles and called the Cayley tree with coordination number 
q and n generations. 

The Bethe lattice may be considered as the interior of 
the Cayley tree. Only sites lying deep inside Cayley tree 
constitute the Bethe lattice. Dealing with the Bethe lat- 
tice we suppose the respective Cayley tree is large enough 
to achieve thermodynamical limit and neglect boundary 
effects. The general structure and topology of the Bethe 
lattice is the same as for Calyley tree. Here we con- 
sider changing along the lattice coordination number q^, 
which is different in general for the different i sites, even 
from the same generation. The Betthe lattice is not uni- 
form now and topological disorder related to it can be 
discussed. 

The Hamiltonian of the Ising model on Bethe lattice 
reads 

H = -j'^SiSj - h^Si, (1) 

where J and h do not depend on the indexes for 
the sake of simplicity. The summation is over all neigh- 
bor Si and Sj spins along the lattice and spin variables 
are Si = ±1. The interpretation of this model in con- 
text of social dynamics is following. The common task 
of social dynamics is the understanding of the transi- 
tion from an initial disordered state to a configuration 
that displays order. The one approach is an analogy to 
physics of the Ising model for ferromagnetic. The Ising 
model may be considered as the simplest description of 
the opinion dynamics. Each spin on the Bethe lattice rep- 
resents a human opinion on some subject (agree/disagree 
units)locating in corresponding lattice site. We describe 
the situation in the simplest sense. The coupling J is pos- 
itive reflecting the fact that the people communicating 
with each other long or even not a very long time usually 
tend to have the same opinion. In other words, humans 
or other social units like communicate with people (so- 
cial units) having usually a good agreement with them. 
However, it is not always true. The discussion of the anti- 
ferromagnetic J is also meaningful. The Bethe lattice in 
our interpretation is not a spatial distribution of humans 
but a construction in relation space and represents pos- 
sible "stable" communications (relations) between "hu- 
mans" Si. Our proposition is to consider "humans" Si 
communicating with (acting on) some finite number of 
neighbors in relation space. This number, which coin- 
cides with the coordination number of lattice qi , may be 
influenced by some external noise and changes in time in 
general. Here we consider qi individually for the each site 
i and the coordination number is not the same for the all 
sites of the same step of the lattice generation. 

An external magnetic field may be treated as an exter- 
nal informational factor, like TV or something similar. 
As we use methods of the statistical physics, the temper- 
ature is introduced and represents the noise in human 
opinion ( measure of uncertainty in taking stable deci- 
sion ) . The main weakness of the model is an absence of 
the small cycles (loops) in relation space. Some effects 



induced by relation loops may be demonstrated on more 
complicated hierarchical lattices (as Husimi lattice). 

III. ANNEALED TOPOLOGICAL DISORDER 

First of all let us recall some general aspects related 
to the spin systems with disorder. Let the Hamiltonian 
may be written as H[J,s\, where J represents disorder 
degree of freedom and s describes the spin variables cor- 
respondingly. Suppose that the spin and disorder degrees 
of freedom are not thermally equilibrated and disorder 
degrees of freedom have their own temperature T. The 
partition function reads: 




DJP{J)exp{-(3F{J)) 



= J DJP{J)exp{^\nZ{J)) 

= J DJP{J){Z{J)r, (2) 

where Z{J) and F{J) are respectively partition function 
and free energy with given disorder J . Here we propose 
to describe the disorder by the probability distribution 
P[J\. We may consider the number n = S as a number 
of replicas within the replica formalism [T5|. The total 
free energy would be 

F = --ln(<Z(J)">d,,), (3) 
n 

where the average over. Using F one would obtain the 
free energy describing quenched disorder after taking the 
limit n — > 0. The n = 1 value corresponds to annealed 
disorder, when both temperatures T , T are equal. The 
value 1 > rt > may be treated here as a "degree of 
quenchness" . The same heuristics arises from the non- 
equilibrium two temperature system viewpoint in [l6i | . 

Now let us consider topological disorder given by the 
set of {qi} variables. We describe only the case of slow 
changing qi compared with the rate of the spin flips. 
This assumption is in agreement with the ferromagnetic 
hypothesis of social opinion dynamics. It is natural to 
suppose that the stable relations dispose to agreement 
between neighbors to a greater extent. So we do not an- 
alyze a more fast rate of qi than it is for the annealed 
disorder. Also we propose the ability of the probability 
factorization P({go, •■•,'7a/}) = P{qa)...P{qM)- 

The one of the best methods developed for the models 
on hierarchical lattices is a so-called dynamical approach 
[l3]-[22|- Here we propose to modify this formalism for 
our purposes. The main quantity to be calculated is 
< Z({5i})" >dis- Let us start from the description of 
the annealed topological disorder n = 1. The main ad- 
vantage of the dynamical approach is that for the models 
formulated on recursive trees the exact recursion rela- 
tions can be derived. When the tree is cut apart at the 
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central site O, it separates into qo identical branches. Ac- 
cording to it for the Hamiltonian ([1]) one can write 



90 



< 9k\qQ,so)...g)^°'{qo,SQ) >dis, (4) 

where independent replicas g'f^ {qo, sq) are invented. The 
Eq. (j4|) also may be rewritten as 



Mo) 



< Z{{q,}) >d^s=J2P'^1o)Y^''P((^'^^o) X 

so 

< 9k{so) >* , (5) 



90 




where 



< 5fc(so) >dis= ^ -P(9o) X! exp(/3JsoSi + /3/isi) 



90 



<g,^^is^)>%;\{6) 



After denoting Xk =< gk{+) > and yk =< 5fe(-) > we 
obtain nonlinear two dimensional mapping 



Xk 



Vk 



/3J+/3h 9-1 I -pj-f3h q-1 



(7) 



where Xk and must be positive. 

As a possible order parameter one may consider mag- 
netization averaged over the disorder. This quantity rep- 
resents averaged opinion in terms of agree/desagree units. 
Using Eq. ^ we get magnetization of the central cite O 
averaged over topological disorder as 



rridi. 



90 



(8) 



Getting fixed points {x,y) of mapping ([7]) and following 
to Eq. ([8|) one would obtain the order parameter rudis 



mdi. 



P{qo) [exp(/?/i)x* - exp(-/3/i)y*] 



P{qo) [eMPh)x''° + exp(-/3/i)y«°] 



(9) 



For the more concreteness we regarded the Bethe lat- 
tice where qi takes values g = 3, 4 with equal probability 
on the each step of the lattice generation and for the 
each site. Corresponding magnetization rridis with re- 
spect to external field f3h at (3 J = 0.5 point is given in 
Fig. ^ . To obtain the magnetization we use the follow- 
ing numerical procedure. At each point (3h, starting from 
Ph — with step A(/3/i) — 0.001, we get fixed points of 
the mapping ([7]). In our simple example, evaluating fixed 
points one have to replace {xk,yk} and {xk-i,yk-i} in 
equations ([7]) on {x,y}. Solution of the obtained system 
of equations with respect to {x,y} gets the set of the 
searched fixed points. In presence of the magnetic field 



FIG. 1: Magnetization curves with /?J = 0.5 in presence of 
topological disorder, where g*, = 3, 4 with equal probabilities. 



there are only one fixed point holding conditions Xk > 
and yk > 0- In case of the zero magnetic field one have 
two stable fixed points and one nonstable, or otherwise, 
depending on temperature. For checking stability of the 
fixed points of ([7]), one can calculate Jacobian of ([7]) as 



J 



e-'+^x- 



2 
2^ 



(10) 

The condition of stability now sounds as | J| < 1. 

It is known that the Ising model on Bethe lattice 
with constant q has a critical point defined by condition 
j3cJ = iln^ [2^. According to [H], it is possible to 
write one dimensional mapping for the partition function 
of the model with fixed q. This mapping has three fixed 
points in zero magnetic field. It is easy to check that the 
temperature corresponding to the case of the three non 
stable fixed points defines by the same condition as the 
critical temperature. When temperature is less then the 
critical temperature, one obtains two stable fixed points 
(spontaneous magnetization) and one non stable x = 1. 
For the temperature greater then critical one, we get two 
nonstable and one stable x = 1 points. This interpre- 
tation of known regimes [1^ around the critical point 
gives a possibility for the numerical investigation in pres- 
ence of the topological disorder. For the evaluation of 
critical temperature of our example g = 3, 4 we propose 
the following numerical procedure. Starting from point 
fiJ = Q and holding = 0, with step A(/3J) = 0.005, we 
check stability of obtained fixed points and looking for 
(3 J at which two stable points become unstable and the 
third one becomes stable. The obtained result is not an 
exact value of the critical temperature, but an interval 
of possible values. The size of this interval depends on 
the precision A(/3J). As a result we found (icJ — 0.725, 
which is greater then the corresponding values for the 
Ising models on Bethe lattice with fixed g = 3 and g = 4. 

Another interesting choice of topological disorder is 
g = 2,3, where g takes corresponding values with prob- 
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ability P{q). Here we consider only uniform probability 
distribution P{2) = P(3) = 0.5. The general study is 
a subject of further research. The coordination num- 
ber q — 2 corresponds to the one dimensional Ising 
model having zero critical temperature. The existence 
and localization of the critical point in this mix of one- 
dimensional Ising model and Ising model on Bethe lattice 
is a subject of special interest. The main idea is an eval- 
uation of the fixed points of corresponding mapping, as 
it was done above for q = 2, 3. However, during this 
procedure one obtains no fixed points holding condition 
x > 0, 2/ > 0. For any positive initial seed of mapping, 
after some finite number of iterations (after 15 — 20 iter- 
ations) Xkjjjk become infinitely big. To deal with these 
big numbers we define new variable Zk = Xk/yk- The 
mapping ([7]) may be represented now as 



Zk = 



P{3){e^'^+'^''zl_^ + e-PJ-^'') + o{n) 



(11) 



The infinitesimals Xk/yf., ^/yk rnay be neglected and we 
obtain that the critical temperature is the same as for 
Ising model on Bethe lattice with constant q — 3. 



IV. QUENCHED TOPOLOGICAL DISORDER 

In case of partly annealed topological disorder < n < 
1 we use the following heuristics. Z quantity may be 
represented as 

z =< zi{q,}r >d^s^ J2 piu^mMph)gll^i+) 



where g^q.y{so) depends on the particular realization of 
disorder. One may obtain recursion relations from the 
previous expression as 



Xk 



Vk 



_ „-(iJ+Ph Ok-l 
- e -^k-l 



l3J-Ph flk-l 



(13) 



denoting Xk = g{qi}{+) and yk = g{qi}{—)- It is assumed 
that on the each step of recursion the power qk takes 
random values with given probability distribution P{q). 
The main idea is to consider that although the mapping 
is stochastic in nature, after some huge number of itera- 
tions one would obtain an attractor set of possible pairs 
{x, y) with corresponding stationary probability distri- 
bution P(x,y). To find the attractor numerically one 
may study the mappings (jl3|) with fixed non stochastic 
values of q^ — q. The fixed points of those mappings 
should be a part of searching attractor. One may start 
stochastic iteration (fTS]) from those fixed points and get 
the set of possible pairs (x, y) . There is a possibility of 
a disconnected attractor (and broken ergodicity), there- 
fore the final set should be based on the summation of 
sets derived from iterations of all possible fixed points. 
During the iteration process one gets not only an at- 
tractor set but also it is possible to evaluate Z. The 
method of iterations takes a chance to calculate value 
of < [eMPh)xl° +exp{-(3h)yl°]'' >d^s■ Although the 
numerical and approximate nature of such evaluation we 
suppose a consistent result in the limit of large number 
of iterations. The magnetization in presence of partly 
annealed disorder may be written as 



jDJP{J)miJ)eM-mJ)) ! DJP{jf-§^{Z{J)r- 



J DJP{J) eM-mJ)) I DJP{J){Z{J)Y 

< [exp(^fe)a;f - exp(-/?/i)yf ] x [exp(/3/i)a;f + exp(-/3/i)y^.°]("-i) >d,. 
< expiPh)xf +eM-WkT >d^s 



(14) 



However, the analysis of mapping (|13p for the simple ex- 
ample where qk = 3, 4 with equal probabilities demon- 
strates some difficulties. After few steps of iteration the 
values of Xk , yk approach to zero or infinity. To over- 
come it we propose a mathematical trick, which should 
give correct values of magnetization for n — s- and exact 
magnetization in presence of quenched topological disor- 
der n = 0. Instead of Xk and yk one may consider new 
variable Zk — Xk/yk and write a one dimensional map- 
ping 

/3J+/3/1 fc-l , -l3J-ph 
Zk = . (15) 



The mathematical trick lies in taking out the means over 
the disorder yk in Eq. (|14p . like it is possible to do with 
constant parameter (real number). It is supposed that 
if yk is infinite quantity or infinitesimal one, than this 
approach is reliable. The magnetization may be rewritten 
now as 
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rudis 



< [exp(/3/i)z^° - exp(-/3/t)] x [exp(/3/i)zf + exp(-/3/i)]("-i) >d. 
< exp(/3/i)zf + exp(-/3/i)]" >dis 



(16) 



r 




0,0 0,2 



The same analysis may be done for the quenched topo- 
logical disorder with q = 2,3. Our study of the sponta- 
neous magnetization shows that the critical point exists 
and the critical temperature differs from zero even when 
the probability of q = 2 is greater than the probabil- 
ity of g = 3 (P(3) = 0.01, 2 < (3cJ < 3). It means 
that the little possibility of branching in one dimensional 
Ising model (rare branches) lead to the critical behavior 
with non zero critical temperature in the thermodynam- 
ical limit. 

As it is easy to see, we mentioned above some proposi- 
tions which are very general (statements about attractor 
and other ) and they need to be studied with mathe- 
matical rigor. However, this statements are useful and 
numerical research gives an evidence of their correctness. 



FIG. 2; Spontaneous magnetization in presence of quenched 
topological disorder, where gj; = 3, 4 with equal probabilities. 



CONCLUSION 



Previous relation Eq. becomes exact in presence of 
quenched topological disorder as it follows directly from 
Eq. (fT4|) . Using this trick is easy to get magnetization 
for our example with g = 3, 4. To obtain magnetization 
we change the values ph like it was done for the an- 
nealed version of disorder. We start evaluation of impor- 
tant means after 10^-th iteration and calculating them 
using the same number of iterations. The result does not 
change if this value is greater, therefore we suppose that 
it is stable. The values of Zk are distributed around the 
fixed points of the mappings ()15p with fixed values of qt ■ 
The dependence of magnetization on external field h in 
presence of quenched topological disorder is given in Fig 

©■ 

To get the position of the critical temperature for our 
example g = 3, 4 we use the basic knowledge about spon- 
taneous magnetization. The dependence of the sponta- 
neous magnetization on temperature provides an oppor- 
tunity to single out the critical point. Iteration process 
and evaluation of the magnetization under the zero mag- 
netic field, as it was described above, give us dependence 
of the spontaneous magnetization on /3J (Fig. The 
critical temperature lies between the critical tempera- 
tures corresponding to the same models on Bethe lattices 
with constant q = 3 and q = A coordination numbers. 



Summarizing written above, we developed an approach 
to deal with topological disorder for the spin models on 
hierarchical lattices. Despite its demonstration only on 
the simple Bethe lattice, a generalization on more compli- 
cated models is straightforward. A wide spectra of prob- 
lems solving within approach of hierarchical models and 
having well physical description within it may be revis- 
ited now in presence of fluctuating coordination number 
of lattice. A complicated disorder may lead to the non- 
trivial solutions of the recursive relations Eq. ([7]). The 
distributions of the Yang-Lee and Fisher zeros in pres- 
ence of disorder is also an interesting problem. Speaking 
about social dynamics we see that the hierarchical mod- 
els take a possibility to simplify the problem within Ising 
paradigm without significant loss of verity and developed 
approach makes it possible to get the major interesting 
quantities. 
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